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Abstract. Howard's conjecture, which states that in the linear instability problem 
of inviscid heterogeneous parallel shear flow growth rate of an arbitrary unstable wave 
must approach zero as the wave length decreases to zero, is established in a mathe- 
matically rigorous fashion for plane parallel heterogeneous shear flows with negligible 
buoyancy force gj3 < 1 (Miles J W, J. Fluid Mech. 10 (1961) 496-508), where j3 is the 
basic heterogeneity distribution function). 

Keywords. Heterogeneous shear flows; linear stability. 



1. Introduction 

In the stability problem of inviscid heterogeneous parallel shear flows, Howard |4|, mak- 
ing use of a novel transformation of the dependent variable was not only able to prove 
the validity of Taylor's conjecture without the restrictive conditions under which Miles 
1 5 1 got the result but he was also able to obtain a semicircular region in the c,-Q-plane in 
which the complex velocity of an unstable wave must lie. Howard also made a conjectural 
assertion to the effect that the growth rate of an arbitrary unstable mode must approach 
zero as the wave number approaches infinity. This conjecture of Howard has also drawn 
the attention of researchers. Banerjee et al | l]| were able to validate the correctness of 
this conjecture for the special case of inviscid homogeneous parallel shear flows. Their 
approach consisted of combining the governing equations and boundary conditions in an 
innovative way and thereby deriving an upper bound of the growth rate under consider- 
ation. Banerjee et al |2| attempted Howard's conjecture in heterogeneous parallel shear 
flows but succeeded only in proving it in the case of the Garcia-type 1 3 1 flows wherein 
the basic velocity distribution has a point of inflexion in the domain of the flow while the 
vertical velocity gradient of the basic density distribution vanishes at this point. 

In this paper we prove, in a mathematically rigorous fashion, that the growth rate of 
an arbitrary unstable wave must approach zero as the wave number approaches infinity 
in the linear instability of non-viscous heterogeneous parallel shear flows with negligible 
buoyancy force (gfi <C 1, |5|) so that second- and higher-order terms in gfi, where j3 is 
the basic heterogeneity distribution function, are negligible as compared to the first-order 
terms in gfi. 



2. Mathematical formulation of the problem 

The basic equations governing the linear instability in a Boussinesq inviscid parallel shear 
flow which is confined between two rigid horizontal boundaries is given by 
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where z £ [21,22] is the real independent variable and stands for the vertical coordinate, 
D = d/dz,U(z) is a twice continuously differentiable function of z and stands for the 
basic velocity distribution, J3(z) is a non-negative continuous function of z and stands 
for the basic heterogeneity distribution, w(z) is the dependent variable and stands for the 
z-component of the perturbation velocity, c = c r + ici is a complex constant in general 
and stands for the complex wave velocity of the perturbation wave with c r as the phase 
velocity and c, as the amplification factor, and a 2 is a positive constant which satisfies 
< a 2 < °° and stands for the square of the wave number of the perturbation wave. 
The boundary conditions associated with the problem are that w(z) vanishes on the rigid 
horizontal boundaries at z = Zi and z = zi i.e., 

w(zi)=w(z 2 )=0. (2.2) 

(The boundaries in the limiting case may recede to ±00.) 

For the existence of a non-trivial solution of eqs (2.1) and (2.2) we have a double 
eigenvalue problem for the determination of c, and c, for the prescribed values of a 2 and 
the flow is unstable if such solutions exist for which the imaginary part c, of c is greater 
than zero. 



3. Mathematical analysis 

Firstly, we prove the following two lemmas: 

Lemma 1. A necessary condition for the existence of a non-trivial solution (w,c, a 2 ) with 
Ci > of eqs (2.1) and (2.2) is that the integral relations 
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are true. 
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Proof. We multiply eq. (2. 1) by w* (the complex conjugate of w) throughout and integrate 
the resulting equation over the domain of z, to get 



J* 2 w*(D 2 - a 2 )wdz - jT 2 w* (^—^j wdz 



I w* ( - 813 ) wdz - 0. (3.3) 
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In order to calculate the first term of the first integral on the left hand side of eq. (3.3), 
we integrate it by parts once and use the boundary condition (2.2), to derive 

jJW + aW)*^ (^--^) MW 0.4, 

Now equating real and imaginary parts of the two sides of eq. (3.4) and cancelling c, (> 0) 
throughout from the imaginary part, we get 

a U"(U-c r ) . |2 
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^{(f/- Cr ) 2 - C 2 } 



k {{U-c r ) 2 + c 2 } 2 
and 

r 2 |2 H o p gP(U-Cr) i ,2 n 

i (£ /-c^ + ^ Hdz - 2 4 { (^-c r)2+C y H = °' 

and this proves the lemma. 



w| z = (3.5) 
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Lemma 2. A necessary condition for the existence of a non-trivial solution (w, c, (X 2 ) with 
Ci > ofeqs (2.1) and (2.2) is that the integral relation 
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are true. 



Proof. We multiply eq. (2.1) by D 2 w* (the complex conjugate of D 2 w) throughout and 
integrate the resulting equation over the domain of z, to get 



jT 2 D 2 w* (D 2 - a 2 )wdz - D 2 w* (jf—^j wdz 

[ Z2 n 2 *( gP 
+ D w 
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wdz = 0. (3.9) 
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In order to calculate the first term of the first integral on the left hand side of eq. (3.9), 
we integrate it by parts once and use the boundary condition (2.2), to derive 



[ Z \\D 2 w \ 2 + a 2 \D w \ 2 )dz- f 1 
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(U-c*) (U-c*) 2 
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On simplification, eq. (3.13) can be written in the form 
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Equating the real parts on both sides of eq. (3.14) we obtain 

r*2 .. „ .„ „. „ „ fz 2 u"(U-c r ) 



(3.10) 



Now from eq. (2.1), we have 

^ = ^+(^)-((^)w. (3.11) 
Taking the complex conjugate of both sides of eq. (3. 1 1), we derive 

»'»'=A- t (^).--( jj A i ).' (3.12, 

Substituting the value of D 2 w* from eq. (3.12) into eq. (3.10), we derive 

U" gj5 



(3.14) 



(3.15) 
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Multiplying eq. (3.1) by a 2 and adding to eq. (3.15) we get 
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which completes the proof of the lemma. 



Theorem 5.1 . If (w, c, a 2 ) is a non-trivial solution ofeqs (2.1) and (2.2) with C{ > and 
heterogeneity factor is small so that second- and higher-order terms in gfi, are negligible 
as compared to the first-order terms in gfi then OlCj — > as a — > °°. 



Proof. Since heterogeneity factor is small so that second- and higher-order terms in gft, 
where j3 is the basic heterogeneity distribution function, are negligible as compared to the 
first-order terms in gfi, therefore eq. (3.8) reduces to 
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Since c, > and \/{{U — c r ) 2 + c 2 } < 1 /c 2 , therefore, from eq. (3.17) we derive that 
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where [(f/") 2 ] max stands for the maximum value of the bracketed expression over the 
interval [zi,Z2]- Further, since (U — c r ) 2 +cj > 2(U — c r )cj, therefore from inequality 
(3.18) we derive that 
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or 



f" 1 {\D 2 w\ 2 +2a 2 \Dw\ 2 )dz 



^.r^jriU)^^, (3 . 20) 

Now, on account of the non-negativity of first integral in (3.20) we get 

cf cf J 

or 

g3c 3 < [(^) 2 Wc,- + ^I^IW. (3 . 22) 
a a 

Inequality (3.22) implies that ac,- -^Oasa^oo. This completes the proof of the theorem. 



Here, it is to be noted that for the case /3 = 0, from inequality (3.22), we get the result 
of Banerjee et al \ 1 1 for the special case of inviscid homogeneous parallel shear flows. 
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